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Abstract 

Wc derived here in a systematic way, and for a large class of scaling 
regimes, asymptotic models for the propagation of internal waves at 
the interface between two layers of immiscible fluids of different densi- 
ties, under the rigid lid assumption and with a flat bottom. The full 
(Euler) model for this situation is reduced to a system of evolution 
equations posed spatially on K d , d = 1,2, which involve two nonlocal 
operators. The different asymptotic models are obtained by expanding 
the nonlocal operators with respect to suitable small parameters that 
depend variously on the amplitude, wave-lengths and depth ratio of 
the two layers. We rigorously derive classical models and also some 
model systems that appear to be new. Furthermore, the consistency 
of these asymptotic systems with the full Euler equations is established. 



Nous ctablissons ici de maniere systcmatique, et pour une grandc 
classe de regimes, des modeles asymptotiques pour la propagation 
d'ondes internes a l'interface de deux couches de fluides immiscibles 
de densite differente, sous l'hypothese de toit rigide et de fond plat. 
Les equations completes pour cette situation (Euler) sont reduites a 
un systeme d'equations d'evolution pose dans le domaine spatial Mr, 
d = 1,2, et qui comprend deux operateurs non locaux. Les divers 
modeles asymptotiques sont obtenus en developpant les operateurs non 
locaux par rapport a des pctits parametres convcnables (dependant de 
l'amplitude, de la longueur d'onde et du rapport dc hauteur des deux 
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couches). Nous ctablissons rigourcuscmcnt des modclcs classiqucs ainsi 
que d'autres qui semblent nouveaux. De plus, on montre la consistance 
de ces systemes asymptotiques avec les equations d'Euler. 

1 Introduction 

1.1 General Setting 

The mathematical theory of waves on the interface between two layers of 
immiscible fluid of different densities has attracted interest because it is 
the simplest idealization for internal wave propagation and because of the 
challenging modeling, mathematical and numerical issues that arise in the 
analysis of this system. The recent survey article of Helfrich and Melville [20] 
provides a rather extensive bibliography and a good overview of the prop- 
erties of steady internal solitary waves in such systems as well as for more 
general density stratifications. The compendium [22J of field observations 
comprised of synthetic aperture radar (SAR) images of large-amplitude in- 
ternal waves in different oceans together with associated physical data shows 
just how varied can be the propagation of internal waves. This variety is 
reflected in the mathematical models for such phenomena. Because of the 
range of scaling regimes that come to the fore in real environments, the 
literature on internal wave models is markedly richer in terms of different 
types of model equations than is the case for surface wave propagation (see, 
e.g. [SI [TO] and the references therein). 

The idealized system that will be the focus of the discussion here, when 
it is at rest, consists of a homogeneous fluid of depth d\ and density p\ 
lying over another homogeneous fluid of depth d<i and density p2 > pi- The 
bottom on which both fluids rest is presumed to be horizontal and featureless 
while the top of fluid 1 is restricted by the rigid lid assumption, which is 
to say, the top is viewed as an impenetrable, bounding surface. This is a 
standard assumption, and is reckoned to be a good one when the pycnocline 
is far from the top, which is when d\ is large relative to the wavelength 
of a disturbance. In the present work, two general classes of waves will 
be countenanced. Both of these require that the deviation of the interface 
be a graph over the flat bottom, so overturning waves are not within the 
purview of our theory (see Figure 1 for a definition sketch). The first, 
which is referred to as the one-dimensional case, are long-crested waves that 
propagate principally along one axis, say along the rr-direction in a standard 
x — y — z Cartesian frame in which z is directed opposite to the direction 
in which gravity acts. Such motions are taken to be sensibly independent 
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of the y-coordinate and can be successfully modeled in the first instance by 
the two-dimensional Euler system involving only the independent variables 
x, z and of course time t. Because the interface is a graph over the bottom, 
these asymptotic models then depend only upon i£l and t, and hence the 
appellation 'one-dimensional'. Among one-dimensional models, the simplest 
are those in which one further assumes that the waves travel only in one 
direction, say in the direction of increasing values of x. Models which we will 
call 'two-dimensional' are not restricted by the long-crested presumption, 
and are consequently more general than the one-dimensional models. They 
are derived from the full three-dimensional Euler system and their dependent 
variables depend upon the spatial variable X = (x, y) 6 M 2 and time t. 

One-dimensional, unidirectional, weakly nonlinear models such as the 
Korteweg-de Vries (KdV) equation , the Intermediate Long Wave (ILW) 
equation [23[ [25] or the Benjamin-Ono equation [5] have been extensively 
used and compared with laboratory experiments [241 [3~Tj [36] . While much of 
our qualitative appreciation of the interaction between the competing effects 
of nonlinearity and dispersion in surface and internal wave propagation has 
been informed by these sorts of equations, they are of somewhat limited 
validity (c.f. [1]). Weakly nonlinear models in two-dimensions have been 
derived by Camassa and Choi [H]. Nguyen and Dias [29] have derived and 
studied a Boussinesq-type system in a weakly nonlinear regime. Fully nonlin- 
ear models were obtained in the one-dimensional case by Matsuno [28], and 
in the two-dimensional case by Camassa and Choi [15] . We mention also the 
interesting paper by Camassa et al. [12] where the aforementioned models 
are compared, in the one-dimensional case, with experimental observations 
and numerical integrations of the full Euler system. In [141 [T5l [28] the anal- 
ysis commences with the full Euler system formulation and the asymptotic 
models are obtained by formally expanding the unknowns with respect to 
a small parameter. It is not easy using this approach to provide a rigorous 
justification of the asymptotic expansion, except perhaps within the setting 
of analytic functions. A different approach has been carried out by Craig, 
Guyenne and Kalisch [T7] in the one-dimensional case. These authors use the 
Hamiltonian formulation of the Euler equations (due originally to Zakharov 
[37] for surface waves and to Benjamin and Bridges [6] for internal waves) 
and expand the Hamiltonian with respect to the relevant small parameters. 
This method provides a hierarchy of Hamiltonian systems which serve as 
approximations of the full Euler equations. Such systems are not always the 
best for modeling, analytical or numerical purposes, however. Indeed, they 
can even be linearly ill-posed in Hadamard's classical sense. In such cases, 
it is necessary in the Hamiltonian framework to proceed one stage further 
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in the expansion, leading to more complicated systems (which may still not 
be well posed). 

The strategy followed here is inspired by that initiated in [8j [9j [TO]. 
Namely, following the procedure introduced in [UJ [TOJ [37] , we rewrite the 
full system as a system of two evolution equations posed on M. d , where d = 1 
or 2 depending upon whether a one- or two-dimensional model is being con- 
templated. The reformulated system, which depends only upon the spatial 
variable on the interface, involves two non-local operators, a Dirichlet-to- 
Neumann operator G[Q, and what we term an "interface operator" H[£], 
defined precisely below. Of course the operator H[£] does not appear in the 
theory of surface waves, and this is an interesting new aspect of the internal 
wave theory. A rigorously justified asymptotic expansion of the non-local 
operators with respect to dimensionless small parameters is then mounted. 
We consider both the "weakly nonlinear" case and the "fully nonlinear" sit- 
uation and cover a variety of scaling regimes. For the considered scaling 
regimes, these expansions then lead to an asymptotic evolution system. As 
in [8j [9j HOj , in each family of asymptotic models may then be inferred 

by using the "BBM trick" and suitable changes of the dependent variables. 
This analysis recovers most of the systems which have been introduced in 
the literature and also some interesting new ones. For instance, in certain 
of the two-dimensional regimes, a non-local operator appears whose analog 
is not present in any of the one-dimensional cases. 

All the systems derived are proved to be consistent with the full Euler 
system. In rough terms, this means that any solution of the latter solves 
any of the asymptotic systems up to a small error. The systems are thus 
seen to be formally equivalent models in terms of the small parameter's that 
arise in the expansions. The advantage of obtaining a family of equivalent 
asymptotic systems is clear from the modeling perspective. One can use the 
flexibility to adjust the linearized dispersion relation to better fit the exact 
dispersion and can choose horizontal velocity variables that are well suited 
to the predictions in view. Mathematically, the choice will be among those 
that are well posed for the particular initial-value or initial-boundary-value 
problem under consideration. When it comes to computer simulation, some 
of the systems are far better suited to the construction of stable, accurate 
numerical schemes and these would naturally be favored. 

The paper is organized as follows. In the next portion of the Introduc- 
tion, the "Zakharov formulation" of the full system is written in dimension- 
less form and the different scaling regimes which will be studied enunciated. 
In Subsection 1.5, a compendium of the outcome of our analysis is offered to 
guide the reader through the rest of the paper. Chapter 2 is devoted to the 
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rigorous asymptotic analysis of the non-local, Dirichlet-to-Neumann opera- 
tor G[C] and the interface operator H[£] mentioned earlier. The asymptotic 
models that result from the use of the expansions of these two operators are 
introduced (and proved to be consistent with the full Euler system) in Chap- 
ter 3. The somewhat technical proof of Proposition 3 is given in Appendix 
A. 

In the present paper, we have refrained from pursuing the analysis to 
the point of obtaining convergence results for the asymptotic systems to 
the full internal waves system. Such a program has been fully achieved in 
the case of surface waves by combining the results of [2] and [ID]. What is 
needed to complete the circle of ideas in the internal wave case is a stability 
analysis of the asymptotic models derived here (that is, an estimation of the 
remainders which comprise the difference between the Euler system and the 
models). Together with consistency, a straightforward analysis would then 
provide a convergence result to the full Euler system, assuming that the 
large time existence results obtained by Alvarez-Samaniego and Lannes in 
[2] for the surface wave system are valid for the internal waves system. The 
latter point is far from obvious; indeed, even the local well-posedness of the 
Euler equations in the two-fluid configuration seems to be an open problem 
in the absence of surface tension (cf [30J for the rigorous derivation of the 
Benjamin-Ono equation for the two-fluid system in the presence of surface 
tension). 

Note finally that the analysis of the present paper could be extended to 
the case of a seabed with structure (a non-flat bottom, see [13] for the case of 
surface waves) and to the case of a two-layer system where the upper surface 
is free rather than restricted by the rigid lid hypothesis (see [U \19\ [28] for 
a derivation of asymptotic models in this situation). These issues are under 
study and an analysis will be reported separately. Of especial interest is a 
comparison of the problem considered with the rigid-lid condition at the top 
and the problem wherein the upper surface is left free in the case where d± 
is relatively large. 
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Notation 

Denote by X the ci-dimensional horizontal variable as described earlier, 
where d = 1,2. Thus, X = x when d = 1 and X = (x, y) when d = 2. We 
continue to use z for the vertical variable. 

The usual symbols V and A connote the gradient and Laplace operator 
in the horizontal variables, whereas Vx. 2 and Ax. z are their d + 1-variable 
version (the gradient in both or all three variables, depending on whether 
d = 1 or 2 and similarly for the Laplacian). For [i > 0, it is very convenient to 
also introduce scaled versions of the gradient and Laplace operators, namely 
V^, 2 = (^V r , d z ) T and A£ f = V Xz • = »A + d*. 

For any tempered distribution u, denote by u or Tu its Fourier transform. 
If / and u are two functions defined on we use the Fourier multiplier 
notation f(D)u which is defined in terms of Fourier transforms, viz. 

f{D)u = fu. 

The projection onto gradient fields in L 2 (M. d ) d is written II and is defined 
by the formula 

vv T 

(Note that II = Id when d = 1.) The operator A = (1 — A) 1 / 2 is equiva- 
lently defined using the Fourier multiplier notation to be A = (1 + jL)) 2 ) 1 / 2 . 
Appearing frequently are the Fourier multipliers T„ and T^ 2 , given by 

T M = tanh^pl) and T M2 = tanh(y7i 2 | J D|); 

where fj,, \ii > 0. 

The standard notation H s (M. d ), or simply H s if the underlying domain is 
clear from the context, is used for the L 2 -based Sobolev spaces; their norm 
is written | • \h s - 

The planar strip S = R d x (-1,0) appears often. The unadorned norm 
|| • || will always be the usual norm of L 2 (S). 

1.2 The Equations 

The Euler system of equations for our system is reviewed here. As in Fig- 
ure 1, the origin of the vertical coordinate z is taken at the rigid top of 
the two-fluid system. Assuming each fluid is incompressible and each flow 
irrotational, there exists velocity potentials <I>j (i = 1, 2) associated to both 
the upper and lower fluid layers which satisfy 

A x , z $i = in nj (1) 



7 



for all time t, where Q\ denotes the region occupied by fluid i at time t, 
i = 1,2. As above, fluid 1 refers to the upper fluid layer whilst fluid 2 is the 
lower layer (see again Figure 1). Assuming that the densities pi, i = 1,2, of 
both fluids are constant, we also have two Bernouilli equations, namely, 

dt^i + hvx^l 2 = --~gz inQi, (2) 
2 Pi 

where g denotes the acceleration of gravity and P the pressure inside the 
fluid. These equations are complemented by two boundary conditions stat- 
ing that the velocity must be horizontal at the two rigid surfaces T± := {z = 
0} and T2 := {z = —d\ — (I2}, which is to say 

d z <S>i = on r h (i = 1,2). (3) 

Finally, as mentioned earlier, it is presumed that the interface is given as 
the graph of a function X) which expresses the deviation of the interface 
from its rest position (A, — d\) at the spatial coordinate A at time t. The 
interface Tt ■= {z = —d\ + A)} between the fluids is taken to be a 
bounding surface, or equivalently it is assumed that no fluid particle crosses 
the interface. This condition, written for fluid i, is classically expressed by 
the relation dtC = yl + | VC\ 2v n-> where v l n denotes the upwards normal 
derivative of the velocity of fluid i at the surface. Since this equation must 
of course be independant of which fluid is being contemplated, it follows 
that the normal component of the velocity is continuous at the interface. 
The two equations 

dtC = \A + |VC| 2 d„$i on r t , (4) 

and 

d n $i = <9 n $ 2 on r t , (5) 

with 

d n :=n-V Xz and n ■= (-VC , 1) T 

v/i + |VC| 2 

follow as a consequence. A final condition is needed on the pressure to close 
this set of equations, namely, 

P is continuous at the interface. (6) 
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1.3 Transformation of the Equations 

In this subsection, a new set of equations is deduced from the internal-wave 
equations ([I])-©. Introduce the trace of the potentials $1 and <3?2 at the 
interface, 

^(i, X) := X, -di + C(t, X)), (t = 1, 2). 

One can evaluate Eq. ((2J) at the interface and use Q and ([5]) to obtain a 
set of equations coupling £ to ^ (i = 1,2), namely 

$C - \A + |vci 2 a„$i = 0, (7) 

4 4 * +,f+ 2 W 2(1 + |VC| 2 ) ) = " P ' (8) 

where in ((7J) and (jSJ), (d n &i) and P are both evaluated at the interface 
= —d\ + Notice that <9 n <3?i is fully determined by tpi since $1 

is uniquely given as the solution of Laplace's equation (fT]) in the upper 
fluid domain, the Neumann condition ([3]) on T\ and the Dirichlet condition 
&i = ipi at the interface. Following the formalism introduced for the study 
of surface water waves in \18\ \19\ [37] , we can therefore define the Dirichlet- 
Neumann operator £?[£]■ by 



g[C]4i = v/i + |vci 2 (a„$i) u= 



-<*!+<' 



Similarly, one remarks that ip2 is determined up to a constant by ip\ since 
<I > 2 is given (up to a constant) by the resolution of the Laplace equation (JT]) 
in the lower fluid domain, with Neumann boundary conditions ([3]) on Y2 and 
d n ^2 = d n <&i at the interface (this latter being provided by ©). It follows 
that ipi fully determines V^2 and we may thus define the operator H[£]- by 

H[C]^i = W>2. 

Using the continuity of the pressure at the interface expressed in ([6]), we 
may equate the left-hand sides of ©1 and (JS])2 using the operators G[C] and 
H[£] just defined. This yields the equation 

W 2 - 7^1) + 5(1 - 7)C + limOH 2 - 7IVV1I 2 ) +AT(CM = 

where 7 = pi /p2 and 

.„> ,s 7(G[C]Vi+VC-V^) 2 - (G[C]^i+VC-H[C]^i) 2 
W M = 2(1 + I VCI 2 ) • 
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Taking the gradient of this equation and using ([7]) then gives the system of 
equations 



d t C - G[C]Vi = o, 

3t(H[Ctyi-7V^i) + </(l-7)VC ( 9 ) 
+ iV(|H[C]Vi| 2 - TlVVil 2 ) + VAT(C^i) = 0, 

for £ and ^l- This is the system of equations that will be used in the next 
sections to derive asymptotic models. 

Remark 1. More precise definitions of the operators G[(] and H[£] will be 
presented in Subsection 11.41 and in Section [2l 

Remark 2. Setting p\ = 0, and thus 7 = 0, in the above equations, one 
recovers the usual surface water-wave equations written in terms of £ and tp 
as in QS1H21EZ]. 

1.4 Non-Dimensionalization of the Equations 

The asymptotic behaviour of © is more transparent when these equations 
are written in dimensionless variables. Denoting by a a typical amplitude of 
the deformation of the interface in question, and by A a typical wavelenth, 
the following dimensionless indendent variables 

A: V . z :- Z ' 



A d\ X/y/gdi, 

are introduced. Likewise, we define the dimensionless unknowns 

r C 7 

a aXyJg/dx 

as well as the dimensionless parameter's 

Pi , di a d\ 

7:=—, 6:=—, e:=—, ji := 

P2 «2 di X z 

Though they are redundant, it is also notationally convenient to introduce 
two other parameter's £2 and fj,2 defined as 

a d\ fi 

Remark 3. The parameters £2 and [12 correspond to e and /1 with c?2 rather 
than d\ taken as the unit of length in the vertical direction. 
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Before writing (|9j) in dimensionless variables, a dimensionless Dirichlet- 
Neumann operator G^[e£]- is needed, associated to the non-dimensionalized 
upper fluid domain 

fii = {(X, z) G R d+ \ -1 + e((X) < z < 0}. 

Throughout the discussion, it will be presumed that this domain remains 
connected, so there is a positive value H\ such that 

1 - e( > H x on R d . (10) 

Definition 1. Let ( G W 2 '°°(R d ) be such that CED is satisfied and let 
-01 G # 3/2 (IR d ). If $i is the unique solution in H 2 (Q.\) of the boundary- 
value problem 



:n) 



^tA$i + fl**i = infii, 
5^i U=0 =0, *i |, = _ 1+eC(JC) = V<i, 

then G^K]V>i e # 1/2 (R d ) is defined by 

G^eQA = -^eVC ■ V*i u= _ 1+e{ + 9 2 <I>i ,_ 1+ef . 

Remark 4. Another way to approach C is to define 

G»[eC]ih = Vl + e 2 |VC| 2 a„$i U= _ 1+e< 

where n $i | z= _ 1+ef stands for the upper conormal derivative associated to 
the elliptic operator /xA$i + <9 2 <I>i. 

In the same vein, one may define a dimensionless operator H M ' 5 [eC]- 
associated to the non-dimensionalized lower fluid domain 

n 2 = {(X, z) G R d+ \ -1 - 1/5 < z < -1 + e((X)}, 

where it is assumed as in (|10p that there is an H 2 > such that 

1 + e5( > H 2 on R d . (12) 

Definition 2. Let ( G I^ 2 '°°(]R d ) be such that {TD} and {TSJ are satisfied, 
and suppose that ip\ G -ff 3 ' 2 (M rf ) is given. If the function $2 is the unique 
solution (up to a constant) of the boundary-value problem 

J ^A$ 2 + 3 2 $ 2 = inft 2 , 

\ d,* 2 \^_ 1/s = 0, d n $ 2 lz= _ 1+si(x) = (1+E 4 c|2)1/2 g^K]^i, 

then the operator H^fe^]- is defined on ^1 by 

H^K]Vi = V($ 2U= _ 1+eC ) G tf 1 / 2 ^). 
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Remark 5. In the statement above, <9 n $2 | z __ 1+eC stands here for the upwards 
conormal derivative associated to the elliptic operator ^A<I>2 + d\ $2> 



vWtt u _ 1+eC = - M eVC • Vd> 2 , z= _ 1+eC + d,$ 2 u= _ 1+s 



The Neumann boundary condition of f)13f) at the interface can also be stated 
as <9 n *2 \ z= _ 1+eC = <9 n $i \ M= _ 1+et - 

Remark 6. Of course, the solvability of (fT3|) requires the condition L d n $2dF ■ 
(where dr = ^/l + e 2 \V(\ 2 dX is the Lebes gue measure on the surface 
r = {z = — 1 + e£}). This is automatically satisfied thanks to the definition 
of G^leQipi- Indeed, applying Green's identity to (fTTI) . one obtains 

[ d n $ 2 dT = [ dn^dT = - [ (/xA$i + d 2 z ^i) = 0. 
Jr Jr Jn L 

Example 1. The operators [£■£]• and H^'*[eC]- have explicit expressions 
when the interface is flat (i.e. when £ = 0). In that case, taking the hori- 
zontal Fourier transform of the Laplace equations ([lip and (|13j) transforms 
them into ordinary differential equations with respect to z which can easily 
be solved to obtain 

GP[0]i/> = -y/Ji\D\ tanh(V7I| J D|)V' and H* 5 [0U = _ tanh (v^l L> l) 

taah(^\D\) 

The equations (jHJ) can therefore be written in dimensionless variables as 

dg--&[eC\^i = 0, 

^(H^^eC]^! - tV^i) + (1 - 7)VC ( 14 ) 
+ |v(|H^K]^i| 2 -7|V^i| 2 )+eV^ 5 (eC,V'i) = 0, 

where Af^ is defined for all pairs (£, -0) smooth enough by the formula 

l(~GntW + VC • V^) 2 - (±G"[Cty + VC • H^[C]^) 2 



2(i + m|vc| 2 ) 



Our work centers around the study of the asymptotics of the non-dimensionalized 
equations (|14p in various physical regimes corresponding to different rela- 
tionships among the dimensionless parameter's e, \i and 5. 

Notation 1. The tildes which indicate the non-dimensional quantities will 
be systematically dropped henceforth. 
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Remark 7. Linearizing the equations (|14p around the rest state, one finds 
the equations 

dtC--G^ 1 = o, 

fi 

9 f (H^[0]^-7VVi) + (l-7)VC = 0. 

The explicit formulas in Example [T] thus allow one to calculate the linearized 
dispersion relation 

, J = {1 } jkj tanh(^I|k|)tanh(^|k|) 

V^tanh^lkl) + 7tanh(^|k|) ' 

corresponding to plane- wave solutions e lk ^~* w *. i n particular, the expected 
instability is found when 7 > 1, corresponding to the case wherein the 
heavier fluid lies over the lighter one. One also checks that the classical 
dispersion relation 

uJ 2 = |k| tanli(- v /7/|k| ) 

y/V 1 

for surface water waves is recovered when 7 = and 5 = 1. 
1.5 Principal Results 

The overall goal here is to propose model systems of equations for the inter- 
nal waves by obtaining the asymptotic form of the equations (|14h in various 
regimes corresponding to different values of the parameters e, 5 and fj,. All 
these asymptotic models are (l+(i)-dimensional systems coupling the surface 
elevation £ to the variable v defined to be 

v :=IF' <5 [£C]V>i -7W>i. (16) 

(For the surface water-wave problem formally recovered by taking 7 = and 
6 = 1, v is the horizontal velocity evaluated at the free surface). We will 
often refer to v as the velocity variable, though its precise interpretation will 
vary. Note that v is essentially the gradient of the second canonical variable 
in the hamiltonian formulation of (|14p . (see for instance [6]). 

It will be rigorously established that the internal-wave equations (|14p 
are consistent with the asymptotic models for (£, v) derived in this paper in 
the following precise sense. 



Definition 3. The internal wave equations (|14p are consistent with a system 
S of d+1 equations for £ and v if for all sufficiently smooth solutions (£, ^1) of 
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(fT4"|) such that §W§ and ([H]) are satisfied, the pair ((, v = H^ 5 [£C]V>i-7'W>i) 
solves S up to a small residual called the precision of the asymptotic model. 

Remark 8. It is worth emphasis that above definition does not require the 
well-posedness of the internal wave equations (|14p . Indeed, these can be 
subject to Kelvin-Helmholtz type instabilities (see for instance [3] and |21j). 
although one might expect a "stability of the instability" result even in the 
face of such instabilities (see |16j). Consistency is only concerned with the 
properties of smooth solutions to the system (which do exist in the classi- 
cal configuration of the Kelvin-Helmholtz problem, even when instabilities 
manifest themselves; see e.g. |33| 132]). In fact, the two-layer water-wave 
system is known to be well-posed in Sobolev spaces in the presence of sur- 
face tension [21]. In consequence, one could simply add a small amount of 
surface tension at the interface between the two homogeneous layers to put 
oneself in a well-posed situation. The resulting analysis would be exactly 
the same and would, in fact, lead to the same asymptotic models. (Such an 
approach is used in [30] for the Benjamin-Ono equation). As the resulting 
model systems do not change, such a regularization has been eschewed here. 

Here is a summary of the different asymptotic regimes investigated in 
this paper. It is convenient to organize the discussion around the parameters 
e and £2 = £<5 (the nonlinearity, or amplitude, parameters for the upper 
and lower fluids, respectively), and in terms of fi and ^2 = (the long- 
wavelength parameters for the upper and lower fluids). Notice that the 
assumptions made about 5 are therefore implicit. 

The interfacial wave is said to be of small amplitude for the upper fluid 
layer (resp. the lower layer) if e <C 1 (resp. £2 <C 1) and the upper (resp. 
lower) layer is said to be shallow if /i <C 1 (resp. fi2 *C 1). This terminology 
is consistent with the usual one for surface water waves (recovered by taking 
px = and 5 = 1). In the discussion below, the notation regime 1 /regime 
2 means that the wave motion is such that the upper layer is in regime 1 
(small amplitude or shallow water) and the lower one is in regime 2. 

1. The small-amplitude/small-amplitude regime: £ ^ 1, £2 <C 1- This 
regime corresponds to interfacial deformations which are small for both 
the upper and lower fluid domains. Various sub-regimes are defined 
by making further assumptions about the size of \i and /i2- 

(a) The Full Dispersion /Full Dispersion (FD/FD) regime: e ~ £2 <C 
1 and /i ~ /i2 = O(l) (and thus 8 ~ 1). In this regime, inves- 
tigated in ^3.1.11 the shallowness parameters are not small for 
either of the fluid domains, and the full dispersive effects must 
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therefore be kept for both regions; the asymptotic model corre- 
sponding to this situation is given in (|26p . 

The Boussinesq / Full dispersion (B/FD) regime: fi ~ e <C 1, 
/X2 ~ 1 . This regime is studied in §3.1.21 and corresponds to the 
case where the flow has a Boussinesq structure in the upper part 
(and thus dispersive effects of the same order as nonlinear effects) , 
but with a shallowness parameter not small in the lower fluid 
domain. This configuration occurs when 5 2 ~ e, that is, when 
the lower region is much larger than the upper one. A further 
analysis of the asymptotic model yields a three-parameter family 
of equivalent systems (see (|2"7|) below). 

The Boussinesq/Boussinesq (B/B) regime: //~//2~£~£2 < Sl- 
In this regime, investigated in §3.1.31 one has 5 ~ 1 and the flow 
has a Boussinesq structure in both the upper and lower fluid 
domains. Here again, a three-parameter family of asymptotic 
systems is obtained (see ([28]) below). 

2. The Shallow Water/Shallow Water (SW/SW) regime: \i ~ // 2 < 1. 
This regime, which allows relatively large interfacial amplitudes (e ~ 
82 = 0(1)), does not belong to the regimes singled out above. The 
structure of the flow is then of shallow water type in both regions; 
in particular, the asymptotic model (see §3.2p is a nonlinear, but non- 
dispersive system, given in (|29p . which degenerates into the usual shal- 
low water equations when 7 = and 5 = 1. It is very interesting in 
this case that a non-local term arises when d = 2. Such a nonlocal 
term does not appear in the one-dimensional case, nor in the two- 
dimensional shallow water equations for surface waves. 

3. The Shallow Water/Small Amplitude (SW/SA) regime: fx < 1 and 
£2 -C 1. In this regime, the upper layer is shallow (but with possibly 
large surface deformations), and the surface deformations are small for 
the lower layer (but it can be deep). Various sub-regimes arise in this 
case also. 

(a) The Shallow Water /Full dispersion (SW/FD) regime: fi ~ e\ -C 
1, e ~ fi2 ~ 1- This regime is investigated in §3.3.11 The disper- 
sive effects are negligible in the upper fluid, but the full dispersive 
effects must be kept in the lower one (see system (j3Tj) below). 

(b) The Intermediate Long Waves (ILW) regime: \i ~ e 2 ~ £2 <C 1, 
fi2 ~ 1- In this regime, the interfacial deformations are also small 
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for the upper fluid (which is not the case in the SW/FD regime). 
This allows some simplifications, as shown in §3.3.21 It is also 
possible (see f|32|) ) to derive a one-parameter family of equivalent 
systems. 

(c) The Benjamin-Ono (BO) regime: n ~ e 2 <C 1, H2 = °o. A 
formal study of this regime is performed in §3.3.31 It is shown 
in particular how to recover the Benjamin-Ono equation as the 
unidirectional limit in the one-dimensional case d = 1. The 
Benjamin-Ono equation is also shown to be a particular case of 
a one-parameter family of regularized Benjamin-Ono equations, 
given in (f3~u) . 

The range of validity of these regimes is summarized in the following 
table. 





e = 0(l) 


£«1 


M = 0(1) 


Full equations 


<5 ~ 1: FD/FD eq'ns 


M< 1 


S ~ 1: SW/SW eq'ns 

S 2 ~ fj, ~ e 2. SW/FD eq'ns 


// ~ e and 5 2 ~ e: B/FD eq'ns 
// ~ e and <5 ~ 1: B/B eq'ns 
S 2 ~ ^ ^ e 2 : ILW eq'ns 
5 = and /x ~ e 2 : BO eq'ns 



Remark 9. The small amplitude/shallow water regime is not investigated 
here. It corresponds to the situation where the upper fluid domain is much 
larger than the lower one, which is more of an atmospheric configuration 
than an oceanographic case. 



2 Asymptotic Expansions of the Operators 

In this section, asymptotic expansions are given of the central operators de- 
fined in the Introduction. The discussion begins with the Dirichlet-Neumann 
operator. 

2.1 Asymptotic expansion of the Dirichlet-Neumann Oper- 
ator G>*[e(]- 

The following lemma connects £ with the vertically integrated horizontal 
velocity via the Dirichlet-Neumann operator G^[eC]-. 
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Lemma 1. Let ( G W 2,OQ (M. d ) be such that [T0\) is satisfied and let ip G 
H 3 / 2 (R d ) and $1 be the solution of (TJ\) with = IfV^ is defined by 



J-l+eC 

then it follows that 

Proof. Let tp G C^°(M. d ) be a test function. Using Green's identity, and with 
the notation of Remark [6j one obtains 



I &[eC[il>ip = [ dn^cpdT 

JR d JT 



Hi 

Defining V^[eQ^j as in the statement of the lemma, it transpires that 



Since the above identity is true for all 99 G C£°(R ), the result follows. □ 

Remark 10. In £ j2.1.1l and ^2. 1.21 below, asymptotic expansions are obtained 
of y /i [eC]V' i n terms of e and (M, respectively. Because of Lemma [H asymp- 
totic expansions of G^[eC]ip then follow immediately. 



2.1.1 Asymptotic Expansion of V*[eC]' when £<1 

When £<1, the approach to obtaining an asymptotic expansion of V^feClV' 
is to make a Taylor expansion in terms of the interface deformation around 
the rest state, viz. 

V»[e()i> = V»[0]iP + e(d (V»[-])()i> + —■ 

(Note, however, that the expansion of V^leQip itself, and not only the con- 
sequent expansion of G^[eQ]ip, is needed so that the elliptic estimate of 
Proposition [3] can be used in the proof of Corollary [1]) . 
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Proposition 1. Let s > d/2 and ( £ H s+3 / 2 (R d ) be such that OH is 
satisfied. Then for tp such that Vip £ H s+1 ' 2 (R. d ), the inequality 

< e 2 C(^-,e^/J2, lCl^+3/2, \Vip\ Hs+ i/2), 

holds for all e G [0,1] and [i > 0, where Tq^ = tanh (^j I l- D l) ^ Ti >At [£] = 
— V / 7o,At(C^o,M^ T )j anc ^ ^ M [ £ C]V' ^ s as defined in Lemma\]\ (so that G^[e(p\p) = 

The key point in the proof is an explicit formula of the derivative of 
the mapping £ i— >• V^[eC]ip, which generalizes the formula obtained in [26] 
for the shape derivative of Dirichlet-Neumann operators. This interesting 
technical point is the subject of the next lemma. 

Lemma 2. Let s > d/2 and suppose that tp is such that Vip £ H s+1 ^ 2 (M d ). 
The mapping H s+3 / 2 (R d ) 3 ( i-» V»[e(]iP £ H s+1 / 2 (R d ) d is differentiable. 
Moreover, for all G H s (M. d ), the derivative of V^fe -]ip at Q in the direc- 
tion C' is given by the formula 

d c (V»[e-}W = -eV»[e(](C'Z»[e(]iP) - eC'(vW - e^VC^KlV), 

where Z»[e(]*P := ^-^{G^e^ + e/iVC • W>)- 

Proof of the Lemma. First, define another Dirichlet-Neumann operator ^[eC]' 
by 

g»H]rP = e z -P»[EC;]V% z % =0 , (17) 

where solves 



(18) 



and where 



(l + e()I dxd e^JI(z + 1)VC 



This operator is the classical Dirichlet-Neumann operator often used for 
the study of the surface water-wave equations and for which an explicit 
expression exists for the derivative of the mapping £ Q^[eQ\tp (see, e.g. 
Theorem 3.20 of [2B] and Theorem 3.1 of [2]). Studying the transformation 
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of the fluid domain into the flat strip — 1 < z < (flattening of the domain) 
reveals that G /i [e£]' ! /' = —G^[—e(]tp (see Proposition 2.7 of [2] and Section 
12.21 below where the same kind of transformation is performed) . It will be 
convenient to consider the operator Q^[— eQ- rather than G^fe^]- because 
this allow us to take over intact some elements of the proof of Theorem 3.20 
in [26] . Moreover, for the sake of clarity, we take e = jjL = 1 in this proof and 
leave to the reader the straighforward modifications for the general case. 
The proof is divided into 5 steps. 

Step 1. One has that G[(]^ = -V • (V[C]), with V[C] = f° x Pi[(]V x,z®dz, 
and where Pi[C] is the d x (d + 1) matrix obtained by taking the last row 
off P[C]- The proof of this result is more or less identical to the proof of 
Lemma [TJ 

Step 2. Denoting by V' the derivative of V[-]tp at £ and in the direction (' , 
one computes 

V = J ^P' I V x ,z$dz + J Pi[C]Vx, g *'dz, 

where Pj and <£' stand, respectively, for the derivative at £ and in the 
direction £' of the mappings £ i— > Pi[C] and £ i— ► <&. 

Step 3. Defining \ = \ z + 1) -j— j- ^ d z &, one has 



° Pi[C]VxA& ~ x) = -v[C](C'-Z[C]^), 
i 



with Z W = To prove this result, first remark that 

w := — x solves the boundary- value problem 

v x,z • p>1 [ £ C}^x,z w = in -Kz<0, 
w u=0 = -C'Z[C]^, 9 z w u= _ x = 0, 



as a consequence of Lemma 3.22 of [26]. The result then follows directly 
from the definition of V[C] - - 
Step 4- The identity 

J (PiVx,z* + Pi[(]Vx,zX)dz = C'(W> - Z[C]il>V£). 

also holds. To establish this, first compute that 

Pfrx** + Pil(]Vx,zX = C'd z {(z + 1)V<&) - VCd z ^ Z+ h 1)2 C'd z ^). 
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The result then follows upon integrating with respect to z. 

Step 5. It now remains simply to put together the pieces. It is deduced from 

Steps 2-4 that 

v = c'(vv - iwvc) - vmc'zicm- 

The result then follows from the observation that if V[(jV is as defined in 
Lemma EEJ then V[(]ip = V[~C]ip- □ 

Proof of the Proposition. A second order Taylor expansion reveals that 

•l 

72 



V»[eC}^ = V»m> + d {V»[e^)C+ / {l-z)di ( -{V^[e-\mX)dz. 

Jo 

Lemma [2] therefore implies that 

V»[£(]ip = V[0]il> 

-eV^0](CG^) - ey^CVV + / (1 - z)d 2 z( -(V^[e-mC, C)dz. 



o 



We saw in Example Q] that G^O]^ = -y^l- ! tanli( v ^|D|)^. Similarly, 
one can check that V^O]^* = tanh |Vj I l p l) The proof of the proposition 
is now clear after appreciating that 

\l-z)d 2 z< (V^[e-mC,C)dz\ H3 <e 2 C(-^,e^l,\(\ Ha+a/2 ,\Vi;\ Ha+1/2 ), 

a fact which is obtained exactly as in Proposition 3.3 of [2]. □ 

2.1.2 Asymptotic Expansion of for Large- Amplitude Waves 

and Shallow Depth (e = 0(1) and ^ < 1) 

For larger amplitude waves, the expansion of the Dirichlet-Neuman operator 
G^leCjip (and also of V^le^tp) around the rest state no longer provides an 
accurate approximation. However, if /i <C 1, which is what we have earlier 
called the shallow water regime for the upper fluid, it is possible to obtain an 
expansion of y /i [eC]V' (and thus of G^eCJV* = \fT^ • ^[^ClVO with respect 
to jU which is uniform with respect to e G [0,1]. 

Proposition 2. Let s > d/2 and ( G J fT s+3/2 (M d ). T/ien /or a// fj, G (0,1) 
and ^ such that V?/> G £P +5//2 (K d ), one has 

Iv^K]^ " M(l " e(W\ HS < fi 2 C(\C\ Ha+3/2 ,\V^\ Ha+5/2 ), 

uniformly with respect to e G [0, 1]), where V^[eQ\tjj is as defined in Lemma 
fj (so that GV[e(]ip = ^//ZV • V^eCW- 
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Remark 11. As in Prop. 3.8 of [2], one can carry out the expansion explicitly 
to second order in fj,, thereby obtaining 

^JIV^sQi; = /i(l - eC)V^ + ^AV^ + 0( M 3 , e^ 2 ). 

Proof. Recall that G^eQip = -Q^-eQij) , where Q^eQ- is defined in <[T7j) . 
and that G^^C]"^ = ~ v^^ 7 " e C]'4' ( see the proof of Lemma[2|). Propo- 
sition 3.8 of [2] shows that 



|^[eC]V - V • (-Ml + eC)VV)L fl < M 2 C(|C|^+3/ 2 , |V^| H 



H5/2; 



An obvious adaptation of the proof shows that the estimate given in the 
statement of the Proposition can be obtained in the same way. □ 

2.2 Asymptotic Expansions of H M ' 5 [eC]- 

Attention is now turned to the interface operator H^' 5 [e^]-. 
The boundary-value problem (|13p plays a key role in the analysis of the 
operator H M ' 5 [e^]-. The analysis of this problem is easier if we first transform 
it into a variable-coefficient, boundary- value problem on the flat strip S : = 
M d x (—1,0) using the diffeomorphism 

s n 2 

° ■ (X,z) i ^ a(X,z):=(X,(l + ed)l + (-l + eC)). 

As shown in Proposition 2.7 of [26] (see also §2.2 of [2]), $2 solves (|13p if 
and only if $ 2 := <3? 2 o~ solves 



with 



V^-^ 2 ^C]V^£ 2 = in 5, 



\ -^£2(2 + ijvc 



(19) 



2 



and where, as before, £2 = sS, fi2 = jz, and V^? 2 = (y^V , d z ) T . 

Remark 12. As always in the present exposition, 9 n $ 2 stands for the up- 
ward conormal derivative associated to the elliptic operator involved in the 
boundary-value problem, 

dn§L2 or = ^ ■ ^ 2 [ £2 C]V^ 2 U =0 or , = -V 
where e z is the upward-pointing unit vector along the vertical axis. 
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An asymptotic expansion of 

H^K^i = V(t 2U=0 ), (20) 

is obtained by finding an approximation §_ app to the solution of (|19p and then 
using the formal relationship H^le^ipi ~ V(§? app \ x=0 )- This procedure is 
justified in the following proposition, whose proof is postponed to Appendix 
lAl so as not to interrupt the flow of the development. The proposition is 
used in both §2.2.11 and §2.2.21 to give explicit asymptotic expansions of 
H^leQipi. To state the result, it is useful to have in place the spaces 

H^(S) = {fGV'(S):\\f\\ Hs , k <^} 

for s £ M and k G N, where = £^ =0 \\ AS ~ jd i /II- 

Proposition 3. Let s > d/2, s > s + 1/2, and ( G H s+3 / 2 (R d ) be such 
that $10\) and are satisfied (the interface does not touch the horizontal 
boundaries). If he H s+1 l 2 > l (S) d+1 and V 6 H s+1 (R d ) d are given, then the 
boundary-value problem 

f V^-Q^[ £2 C]V^ = V^-h »5, (2i) 
\ d n u\ z ^ Q = ■ V + e z • h u=0 , d n u\ z= _ t = e z • h u= _ 1 

admits a unique solution u. Moreover, the solution u obeys the inequality 

|v« u=0 | H . < -^c(^, e r x ^rMci^ + 3/ 2 )(iih|i HS+1/2 , 1 + 

uniformly with respect to e 2 G [0, e™^] A*2 G (0, n™ ax ). 

Remark 13. In the case of a flat interface (£ = 0), Example Q] shows that 
±G^[0]V>i = • V with V = Xtanh^lDl)^!. Consequently, (fT9|) . 

(I20p and Proposition [3] (with h = 0) show that 



s + 1 . 



which is exactly the estimate one could have deduced from the explicit ex- 
pression for H^[0]- given in Example [1] (except that using the latter ap- 
proach gives an estimate in H s rather than in H s+l . The £P-type result 
does not in fact carry over to the general case of non-flat interfaces). 

Remark 14. Suppose we take h = and V = V^fe^]^ m Proposition 3. 
By Lemma [U one has G^leQijj = ^/JIV ■ V^le^ip, and so it follows that 
Viii =0 = H /i ' <5 [e£]?/>. The Proposition thus provides an estimate of the 
operator norm of H^fe^]. 
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2.2.1 The Small- Amplitude/Small- Amplitude Regime: e <C 1, 

e 2 « 1 

In this regime, it is assumed that the interface deformations are of small 
amplitude for both the upper and lower fluids. The asymptotic expansion 
of the operator H^^leQ is thus made in terms of e and e 2 = We proceed 
by first constructing formally an approximate solution 3>, to (fl"9j) in the 
form 

ta PP = $ (0) +£2$ (1) . 

This formal approximation is then justified rigorously in Corollary [1] below. 
We may write from the expression for Q^^CL 

■ M Vg, = Ag, + e 2 V^ • QiVjJ?, + eiVg, • Q 2 V^, 

with 

Ql ~ { -jjr 2 {z + i)vc T -C 

and 

/ \ 

<32 = n C 2 +M2(^+1) 2 |VC| 2 
1+E2C 



It follows that 

+ e 2 (A^ • $W + • QiV^ z $(°)) + 0(e|). 

Similarly, we obtain 

ft^twp = ^ (0) +£2(e. -QiVg^W 1^ +0{el). 

Since it is known from Proposition [1] that 

- 5 &[eC]^i = ■ (T , M V^i) + e 2/ u 2 V • (-C + T 1>M [C])V^i + 0(ie^ 2 ), 

one therefore deduces that 3> app solves (fT9j) up to order 0(^2 + ji^\^2) 
provided that <J>(°) and $W solve 

r A£ 2 z $(°) = 0, 

1 u=0 = • (r , M v^i), d z m u= _ x = 0, 
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which is obviously solved by d>(o) (X ,z) = 
and 

\ d z $V lz=0 =A, 8,^^=0, 

with A = fj, 2 V ■ (-C + TtJCWi ~ e z • QiV£%$(°) u=0 . Because -V£% • 
QiV^$(°) = A^J(z + 1)C^$(°)] and 

A = /i 2 V • [ - C + -TUClWi] + M2 V • (CV^W) + + l)C9x$ (0) ) | l=0 , 

it results that = (z + l)<^d z ^ ^ + u, where u solves the boundary-value 
problem 

A X,z u = °» 

8 z u ]z=0 = / u 2 V- [-C + T liM [C]V^i] +^ 2 V-(CV<I>( )) 5^1^ = 0. 

This latter boundary- value problem can be explicitly solved by taking the 
Fourier transform in the horizontal variables and solving the resulting ordi- 
nary differential equation in the variable z. One obtains from this calculation 
that 

V7 n- !gj rrr^i i tanh(^I|D|) 

Vn l- = ^ tanh(^|^|) n[C(1 + tanh(^|^|) )V ^ l] 
^ r tanh^lgp ^ tanh(^|D|) 

+ ^v[ tanh(v _ |Z)|) (C m ] ' 

Since $W = n + ( 2 + l)(d z <f>^ and 

v[(, + i)c^°']i„ =VMiv(c tMh ;ff |c|) A^), 

it is deduced immediately that V^ 1 ) i = B((,Vipi), where 

+ ^[(l + ^^)(C^^A*)].(22) 

The rigorous result concerning the asymptotic expansion of the operator 
H^ ,<5 [e£] in the present regime, which is a corollary of Proposition [3l may 
now be stated and proved. 
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Corollary 1 (Full dispersion/Full dispersion regime). Let to > d/2, s > 
t + 1/2, and ( G H s+3 / 2 (R d ) be such that J22]) and J2IP are satisfied. Then, 
for all ijj x such that V^i G £P +5/2 (M d ), 



7i s 



c^tc , j- rraai „mai „moi iai MTV,,/, I 



where the bilinear mapping B(-, •) is defined in This estimate is uni- 

form with respect to e G [0,1], ^ G (0, /i max ) and 5 G (0,<5 max ) smc/i that 

,, A* c (C\ ,,max\ 

Proof. The computations above show that 

v^-o^Mv^t^ = e |v^ -h, 

with h = QiVx,z® {1) + ( 32V^ Z ($(°) + e 2 $ (1) ). It is also easy to check that 

5ntap P | z=0 = V^v-(r 0iM v^ 1 )+e 2/ u 2 v-(-c + r liAt [c])v^i + ^e z -h u=0 , 

^nilapp U=-i = £2 e z- h U=_i- 

Therefore, the difference v = §L app — 5j? 2 satisfies the boundary- value problem 
V% z -Q^[e2QV% z v = elV% z h, 

<9„u | z=0 = y/JpiV ■ V + e|e, • h u=0 , d n v \ z= _ t = e|e z • h| z= _ l5 

with V = (T 0iM V^i) + e 2v ^2"(-C + T liM [C])V^i - V^K]V>1, and where 
V^eC]^! is given by Lemma [TJ Applying Proposition [3] in this situation, it 
is immediately deduced that |Vu| z=0 |#s is bounded from above by 

c'(^,^ j/i r B jci^3/ a )(-|=iihii^i/ a ,i+-^iF|^i). 

H 2 v/M2 J 1*2 



The stated result is thus a direct consequence of Proposition Q] and the 
observation that ||h|| H .+i /3 ,i < C(^, <5 max , /x^, |C| HS+3 / 2 )| V</>i | ffs +3/2. □ 

This section concludes with two specializations of Corollary [1] that obtain 
when additional smallness assumptions are made on the parameters fi, /U 2 
or on 5. These simple consequences of Corollary Q] will be useful presently. 
The two additional regimes we have in mind are the following. 
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1. The Boussinesq/Full dispersion regime. This regime is obtained by 
assuming that fi ~ e and /i 2 ~ 1 (and thus 5 ~ e 1 / 2 ) in addition to 
the assumptions e<1 and £2 1 which are required if one wants 
Corollary [I] to provide a good approximation. 

2. The Boussinesq/Boussinesq regime. Here, it is assumed in addition to 
e <C 1 and £2 <C 1 that fj, ~ e and /i2 ~ £2 (and thus 5 ~ 1). 

Corollary 2 (Boussinesq/Full dispersion regime). Let tQ > d/2, s > to + 
1/2, and ( G H s+3 / 2 (R d ) be such that (Q2|) and (OJP are satisfied. Then, for 
all tpi such that VV>i G iP +5 / 2 (R d ), £/ie inequality 

|H^[eC]^i - vW coth(V^2~|r>|) [ - VVi - f AV</>! + eII«V^i)] |^ 

p2 _i_ _2 11 
^ / 2 ' c 1 _, 1 _ , ,1/2 r\ /~t( rmax ,,max ,.max iai MV7„/, I 

v/^2 -Hi -"2 



where II = — tjjff > holds uniformly with respect to e G [0, 1], fi £ (0,/i max ) 
and 5 G (0,o" maa: ) suc/i that fi 2 = f £ (O,/^)- 

Remark 15. When e -C 1, /i ~ s, [i% ~ 1 (and thus b ~ e 1 ^ 2 ), the three 
components of the error estimate are all of the same size 0(e 2 ). 

Proof. The result is obtained by using tanh(- v //^l^'l) ~ — fi^/JI^\D\ 3 

when [i is small in Corollary [TJ □ 

Similarly, one may also deduce from Corollary [1] the following result in the 
Boussinesq-Boussinesq regime. 

Corollary 3 (Boussinesq/Boussinesq regime). Let tQ > d/2, s > to + 1/2, 
and C G H s+3 / 2 (R d ) be such that (EJ) and are satisfied. Then, for all 
ipl such that Vipi G H s+5 ^ 2 {M. d ), we have 

6 M _ 1 

< $±£. + f + £ 2 )C(-^, -L -i-, |C|^ + 3/ 2 )|V^i|^ +5/2 , 

where IT = — . Moreover, this estimate is uniform with respect to £ £ 
[0, 1], fi£ (0,/i max ) and 5 G (5 min ,5 max ). 

Remark 16. When £~£2~/i~/U2^Cl (and thus 5 ~ 1), the last two 
components of the error estimate are of size 0(e 2 ), but the first is of size 
0(e 3 / 2 ). This loss of precision is not seen at the formal level. It comes from 
the l/^yji2 term in the elliptic estimate provided by Proposition El 



|H^[<]^ - ( - «5VVi - -/i(l - t^AV^ + e 2 (l + <S)II(CVVi)) L 
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2.2.2 The Shallow- Water/Shallow-Water Regime: /i 2 «l 

In this regime, large amplitude waves are allowed for the upper fluid (e = 
O(l)) and for the lower fluid (£2 = 0(1)). Assuming that fx <C 1 and 
fj>2 *C 1 raises the prospect of making asymptotic expansions of shallow- 
water type, in terms of \i and [12- As before, the plan is to formally construct 
an approximate solution & app to (fT9l) having the form 

The formal approximation is then rigorously justified (Corollary [5] below) 
and the desired expansion results. From the expression for Q^ 2 [e^CL we may 
write 

n-2 

with \i2 = 1 + £2C an d 

h 2 Idxd -e 2 (z + l)V(\ 

S2(Z + 1)VC T £i(Z+ g' VC ' 2 J' 

It follows readily that 

Similarly, one infers that at z = and z = — 1, 

= t^ (0) + M2(e z • QiVx,,<I> {0) + ^-d-M l) ) + O^l). 
ri2 1^2 

Since it is known from Proposition [2] that 

1 u 2 

-G^eC^i = W ■ (/11VV1) + O(^) 



(with hi = 1 — eC), it is clearly the case that solves (|T9|) up to order 
0{\i\ + provided that $( 0) and $W solve 
f 52^,(0) = Qj 

1 a^w u=0 '= 0, a,*w> = 0, 

which is obviously solved by any &(°\X,z) = $>(°\X) independent of z, 
and 

f = -h 2 2 A&°\ 

{ d z <s>w u=o = ^ 2 ( £2 vc • v$(°) + av • (/ti v^i)) , d z $w u=-i = 0, 
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where we have used the fact that does not depend on z. Solving this 
second order ordinary differential equation in the variable z with the bound- 
ary condition at z = yields (up to a function independent of z which we 
take equal to for the sake of simplicity) , 

^ (1) = -Y^Ad>(°) + z(a^ lu=0 ). 

Matching the boundary condition at z = — 1 leads to the restriction 

V • (/i 2 V$ {0) ) = -SV • (hiVipi), 

which implies that II^V^ )) = II(— ShiVi/ji), where II = — T^js is the or- 
thogonal projector onto the gradient vector fields of L 2 (R d ) d defined earlier. 
We will solve this equation thanks to the following lemma. 

Lemma 3. Assume that £ G L OD (R d ) is such that |e 2 C|oo < 1- Let also 
W G L 2 (R d ) d . Then 

i. One can define the mapping 0[e 2 Q as 

L 2 (R d ) d L 2 (R d ) d 

Q[e2C] : u i — > ^(-i) n (n(e 2 cn-)) n (n^7) 

n=0 

ii. There exists a unique solution V G L 2 (R d ) d to the equation 

V • (h 2 V) = V • W, (h 2 = l + e 2 Q 
such that UV = V and one has V = Q.\_E2(,\W ; 

Hi- If moreover Q G H s (R d ) andW G H s (R d ) d (s > d/2+1) then 0.[e 2 C,]W G 
H s (R d ) d and 

\Q[e2<]W\ H ' < C(\e 2 C\h;-i ^^)\W\h^. 

i - l^2<,|oo 

Remark 17. In dimension d = 1, one has II = 1 and the first point of the 
lemma simplifies into V = j^W so that the proof is trivial. 

Proof, i. The result follows from the observation that under the assumptions 
of the lemma, one has 

linden-) \\ L ^ L 2 < |e 2 c|oo < i, (23) 
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so that the series used to define Q[e2(]U converges in L 2 (R d ) d . 

ii. Let us first check that V = £l[£2(]W is indeed a solution of the equation 
stated in the lemma. Since V = IIV, one can remark that 

V-(e 2 (V) = V-(n(e 2 CIBO) 

oo 

= -v-^(-ir(n( £2 cn-)r(nw) 

n=l 

= -V ■ (V -ILW), 

from which one deduces easily that V • (I12V) = V • W. 
Let us now turn to prove uniqueness of the solution by proving that one has 
necessarily V = if VF=0. To check that this is the case, just remark that 
from the equation V • (I12V) = and the requirement that HV = V, one has 

v = -n(e 2 cm/) ; 

since ||n(e2Cn^)||L2^ i 2 < |£2C|oo < 1, it follows easily that V = 0. 

iii. It is clear from (f23j) that IQ^Cl^k < 1 _| £ 1 2 ^| ■ Now, applying A s to 
the equations, one gets 

V • (h 2 A s V) = V-W, 

with W = A S W + [A S ,£2C]^- The result follows therefore from the L 2 
estimate, a standard commutator estimate and a simple induction. □ 

If Lemma [3] is applied with V = V$' ', W = — <5/iiV - 0i, there results 
the equation 

V$ {0) = -5 a[e 2 C]OiVV>i)- 
Note that when d = 1, this reduces to 

n,2 

The following corollary of Proposition [31 which gives the needed asymptotic 
expansion of the operator EP^fe^] in the present regime, now comes into 
view. 

Corollary 4 (Shallow water/Shallow water regime). Let to > d/2, s > 
t + 1/2, and ( £ H s+3 / 2 (R d ) be such that ^ and HE) are satisfied. Let 
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hi = 1 - e( and h 2 = 1 + e 2 Q and let fa be such that Vfa G H s+5 / 2 (M. d ). 
Then it follows that 

\H^ 6 [e(}fa + S Q[e 2 (](hiVfa)\ H , 

< ^+/i 2 )c((i-,5(i-F 1 ))- i ,^, e r x ^rMci^ + 3/ 2 )ivvii^ +5 / 2 , 

ti 2 

uniformly with respect to e G [0,1], fi G (0,1) and (5 < j^j^ such that 
e 2 = e5 G [0,e^ ax ] and ^ 2 = £ G (O,//^). 

Remark 18. When e~£2 rs -'A i ~A i 2^1 (and thus <5 ~ 1), one deduces from 
the above corollary that U^' S [eC]ipi = —SVfa + 0(e), which is consistent 
with the asymptotic expansion provided by Corollary [3J A similar matching 
would have been observed for the next order terms if we had computed them 
in Corollary |H 

Remark 19. When d = 1, one has 5 £2[e 2 C](^iV^i) = 5j^d x fa. 

Proof. Since (fTUD , CE2D and the condition 5(1 -Hi) < 1 imply that |e 2 C|oo < 
1, one can use Lemma [3] and the computations above indicate that 



with 



ff 2 v£, ■ h, 



It is also easy to check that 

dntapp |* =0 = • (hiVfa) + fi 3 2 /2 e z ■ h u=0 , 

3/2 

9n!kappU=-i = N e z h U=-i- 

Thus, the difference u = 5^ — <£ 2 satisfies the boundary- value problem 
f V^-Q^hC]V^u = ^ /2 V^ -h, 

\ d n u\ z=0 = iffP2~V ■ V + ^2 /2 e z • h| a=0 , 5 n -U| z= _ 1 = n 2 /2 e z • h| z= _ r 

with V = /ii<5 \/Ji 2 V 'fa — V^leQfa, where y^ l [eC]V'i is given by Lemma [TJ 
One concludes from Proposition [3] that |Vui =0 \h 3 is bounded from above 
by the quantity 
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The result is a direct consequence of Proposition [2] since 

1 1 h H^+l/2,1 < SC( — , E™ aX , fl™ aX , \(\ H s+3/2)\Vlf)l\ Hs +3/2 
n-2 



□ 



2.2.3 The Shallow- Water/Small- Amplitude Regime: /i<l, £2<1 

It is now presumed that both fi and £2 are small, but no such restriction is 
laid upon e nor /J2- So, this regime is not a subcase of the regimes inves- 
tigated in Sections 12.2.11 and 12.2.21 We construct an approximate solution 
Q-app t° (|19p exactly as in §2.2.11 but only a first-order approximation of the 
form 

$ =$(0)- 

—app ! 

will be required. Since /x 1 here, Proposition Q] may be utilized to write 
]G*[eCWi = fr-(hiVih)+0(!j). 

Just as in £ j2,2.H it can be shown that must solve the boundary- value 
problem 

\ d z $w u=0 = f v • (/nWi), u=-1 = 0, 

which is to say that 

(0) cosh( v ^ 2 "(2: + 1)\D\) 1 _ /, _ M 

cosh(y7l2|Z)|) taoh.(y7^|Z)|) y 

The following result is proved exactly as was Corollary [TJ 

Corollary 5 (Shallow water/Small amplitude regime). Let io > d/2, s > 
t + 1/2, and C e i? s+3 / 2 (M d ) 6e sucft i/W |22]) arid are satisfied. Then, 
for all V>i swc/i that Vipi G i? s+5 / 2 (]R <:( ), it is tae case taat 



|h^ 5 K]^i + Vv\d\ coth(v^|D|)n(/ ll Wi) 



/J + £2y/Ji . 1 1 j-max ,,max \/-\ MV7„/, I 

V^2 -"I ^2 



where h\ = 1 — e£ and II = — is given by \22$. This estimate is 
uniform with respect to e G [0,1], /i G (0,1) and 5 G (0, £ ma:r ) swc/i t/iat 
— ^ G (U, H 2 ). 
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Remark 20. Several regimes fall within the range of Corollary 

• The SW/FD regime: when \i <C 1, £2 *C 1 and e ~ /i 2 ~ 1 (and thus 
5 2 ~ ~ el); the precision of the approximation is O(fi). 

• The ILW regime: if /x ~ e 2 <C 1 and /li 2 ~ 1 (and thus 5 2 ~ /z ~ e 2 ); 
in this case, the estimate in the corollary can be simplified without 
adverse effects on the precision of the approximation to simply 

H^[eC]Vi = -yfc\D\ t5oth(v^i|D|)V^i + 0(fi). (24) 

• The BO regime: if /i <C 1 and 5 = (and thus /x 2 = 00, £2 = 0), one 
gets formally from (124p that 

lP' 5 [ e C]V>i ~ -VmPIV^i. (25) 
3 Asymptotic Models for Internal Waves 

The preliminary analysis in Section 2 allows us to derive the various asymp- 
totic models referred to in the Introduction. 

3.1 The small amplitude/small amplitude regime: e < 1, 

e 2 < 1 

Derived first are various models corresponding to the case wherein the in- 
terface deformation are small for both fluids. Different systems of equations 
obtain, depending on the sizes of the paramaters e, \x and 5 (and thus £2 
and fj, 2 )- 

3.1.1 The Full Dispersion/Full Dispersion Regime: e ~ e 2 1 

and fj, ~ yu 2 = 0(1) 

An asymptotic model can be derived from (|14p by replacing the operators 
and FP ,<5 [e£] by their asymptotic expansions, provided by Proposition 
[T]and Corollary Q] in the present regime. The following theorem shows that in 
the present regime, the internal wave equations are consistent with following 



32 



FD/FD system; 

I _£2_ V . / T M T M2 P(V T M2 

+ v^M ^T (l2 +T, D i'.' 7 T, ;! +T/ii 

" eV • (C^r;v) + e\D\T,(C^ • (^v)) = 

dtv + (i - 7 )vc 

(26) 

where as before, T M = tanh(y / /7|Z)|), T M2 = tanh(y//l2|.D|) an d the bilinear 
mapping B(-,-) is given by (f22|) . 

Theorem 1. Let < 5 mm < 5 max . The internal waves equations p4\ ) are 
consistent with the FD/FD equations 1126}) in the sense of Definition^ with 
a precision 0{e 2 ), and uniformly with respect to e € [0, 1], fj, € (0, n max ) and 



Remark 21. One can give a more precise estimate of the precision, as in 
Corollary Q] for instance. It simplifies the exposition to use the notation 
0(e 2 ) and the associated rough estimate of the precision. We follow this 
policy throughout the discussion. 

Remark 22. It is straightforward to check that the dispersion relation of 
(|26p is exactly the same as (|15p . which is the reason we refer to (|26p as a 
"full dispersion" model. In particular, (f26j) is linearly well-posed provided 
that 7 < 1. 

Proof. First, notice that with the range of parameters considered in the 
theorem, one has e ~ £2 when e — > 0, while fi ~ = 0(1)- By the 
definition (|16p of v and using Proposition Q] and Corollary [IJ one deduces 
from ([HI) that 



d t( ~ 7PW\ • ( T m V ^) + eV • (Wi) " eP|T M (C^| • (T M W0) = 0(e 2 ) 
<9 4 v + (1 - 7 )vc 

+ |V(|H^K]^i| 2 - 7lW>i| 2 ) + e^V(^ • (T^V^i)) 2 = 0(e 2 ). 
It follows from Corollary [1] and the relation H^' 6 [eQtpi = v + jVipi that 

^ = ~ ^fr- ( v + £2B(C ' vr^Pr- v) ) + ° (e2) - 

The result is view is now apparent. □ 
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3.1.2 The Boussinesq/Full Dispersion Regime \i ~ e <S 1, fj,2 ~ 1 

We show here that in this regime (for which one also has 5 2 ~ e and thus 
e 2 e 3 / 2 1) ; the internal waves equations (|14p are consistent with the 
three-parameter family of Boussinesq/FD systems 

(l- M 6A)d t C+-V-((l-<M 

- ^\D\ coth(^|D|)V • V/9 + ^ (a - coth 2 (^|D|)) AV • V/3 = 

(1 - fjLdA)dtvp + (1 - 7 )VC - 7fv| V/3 | 2 + /xc(l - 7 )AVC = 0, 

(27) 

where vg = (1 — fi(3A) x v and the constants a, 6, c and (i are defined now. 
Theorem 2. Let < c min < c max , < < fxf ax , and set 

a = -(1 — ai — 3/3), b = —ai, c = (3ct2, d = (3(1 — a^), 
3 3 

with OLi > 0, (3 > and 02 < 1. these choices of parameters, the in- 

ternal wave equations |J^| ) are consistent with the Boussinesq/FD equations 
(21) in the sense of Definition^ with a precision 0(e 3 ^ 2 ), and uniformly 



with respect to e S [0,1], (J, G (0, 1) and 5 G (0, 1) satisfying the conditions 

c min < i < c max ^ mm < M < max_ 

- fi ~ ^ ~ 6 2 ~ 1 

Remark 23. The dispersion relation associated to (|27p is 



2 1-7, l2/ „ l2 l-^|k|coth(^|k|)-/i|k| 2 (a-4coth 2 (^|k|)) 

w = ~T~ |k| 1 } (i + ^|k| 2) (i + Hk| 2 ) 

and (|27p is therefore linearly well-posed when b, d > and a, c < 0. Notice 
that in the case a\ = 012 = /3 = 0, one has = 3 and b = c = d = and 
the corresponding system is thus linearly ill-posed. The freedom to choose a 
well-posed model is just one of the advantages of a three-parameter family 
of formally equivalent systems. The same remark has already been made 
about the Boussinesq systems for wave propagation in the case of surface 
gravity waves [HI flU]). 

Proof. The proof is made in several steps, corresponding to particular as- 
sumptions about the parameters 01,02 and f3. Throughout, use will be 
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freely made of the relations /i ~ e and fj,2 ~ 1. 



Step 1. The case a\ = 0, j3 = 0, a.2 = 0. From the expansion of the 
Dirichlet-Neumann operator, see Remark [TTT it follows as in the previous 
section that 

d t C-V- ((1 - eC)VVi ) - f V ■ AV^ = 0(e 2 ) 
a t v + (1 - 7 )VC + |V(|H^[eC]^i| 2 - 7lW>i| 2 ) = 0(e 2 ), 

where the fact that 0(/jl) = 0(e) has been used. From the relation H At,<5 [££]?/>i 
v + ryVipi , and Corollary [21 it is seen that 

vv* = - -v - [1 + £ A - £ n(c-)] vy>i + o(e 2 ). 

7 7 J -/i2 

Again using the fact that 0(/x) = O(e), one concludes that 

„ , 1 JJi \D\ a A 3. 

7 7 2 T M2 7 J T£ 2 

and the result follows. 

Step 2. The case ai > 0, = 0, ai = 0. We use here the the classical 
BBM trick [7] . It is clear from the first equation that 

a i c = --v-v + o(e 1 / 2 ), 

7 

from which it is inferred that 

V • v = (1 - ai)V • v - a^d t C + 0(e 1/2 ). 

Replacing V • v by this expression in the component g^AV • v of the first 
equation of the system derived in Step 1, leads to the desired result. 
Step 3. The case a.\ > 0, (3 > 0, a<i = 0. Replacing v by (1 — /x/5A)v / g in 
the system of equations derived in Step 2, and neglecting the 0(e 3 / 2 ) terms 
is all that is required in this case. 

Step 4. The case a± > 0, > 0, «2 < 1- We use once again the BBM trick. 
From the second equation in the system derived in Step 3, one obtains that 
for all Q!2 < 1) 

dtvp = (1 - a 2 )d t V(3 - a 2 (l - j)V( + 0(e). 

If this relationship is substituted into the system derived in Step 3, the result 
follows. □ 
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3.1.3 The Boussinesq/Boussinesq Regime £~ / u~£2~/^2^1 

In this regime, the nonlinear and dispersive effects are of the same size for 
both fluids; the systems of equations that are derived from the internal waves 
equations (|14|) in this situation are the following three-parameter family of 
Boussinesq/Boussinesq systems, viz. 

(l - //6A) d t C + ^-^V • + e i-=^V • (Cvp) + iiaV • A V/3 = 

< r 2 

(l - fMiAj d t v p + (1 - 7 )VC + | ((5+ ~ 7 ] 2 V|v^| 2 + ^cAVC = 0, 

(28) 

where va = (1 — fx/3A) x v, and where the coefficients a, 6, c, d are provided 
in the statement of the next theorem. 



1+7^ 



_ (l-ai)(l+7<*)-3fl?(7 + <5) 
a " 35(7 + 5)2 ' 1 3*(t+«) ! 

c = (3a 2 , d = /?(l — a 2 ), 

TOi/i ai > 0, (3 > and a 2 < 1. Witft £/ws specification of the param- 
eters, The internal wave equations (7^1) are consistent with the Boussi- 



nesq/Boussinesq equations Ii28\) in the sense of Definition^ with a precision 
0(e 2 ), and uniformly with respect to e £ [0,1], fj, € (0, 1) and 5 £ [5 min , S max ] 
such that c min < % < c max . 

Remark 24. Taking 7 = and 5 = 1 in the Boussinesq/Boussinesq equations 
(f28l) . reduces them to the system 

(1 - ^A)d t C + V • ((1 + eC)v) + ix — ^p^ V • Av = 
(l - n/3(l - a 2 )A)5 t v + VC + |v|v| 2 + M/3a 2 AVC = 0, 

which is exactly the family of formally equivalent Boussinesq systems derived 

m [sung. 

Remark 25. The dispersion relation associated to (|28|) is 

a (^-Hk| 2 )(l-7~Hk| 2 ) 

W 1 1 (l+/i6|k| 2 )(l + Hk| 2 ) ' 

It follows that (|28p is linearly well-posed when a, c < and b,d > 0. The 
system corresponding to ct\ = a 2 = (3 = is ill-posed (one can check 



3G 



that a = 3g?~^p > 0). This system corresponds to a Hamiltonian system 
derived in [17] (see their formula (5.10)). As mentioned before, the present, 
three-parameter family of systems allows one to circumvent the problem of 
ill-posedness without the need of taking into account higher-order terms in 
the expansion, as in [IT]). 

Proof. The proof is again made based on various possibilities for the param- 
eters in the problem. For this regime, we have that e ~ fj, ~ £2 ~ ^2 as 
e — > 0. The overall idea of the argument is the same as evinced in the proof 
of Theorem [TJ 

Step 1. The case ot\ = 0, /3 = 0, oli = 0. Using Remark II II and Corollary [3] 
(instead of Proposition [1] and Corollary [T] as in the last theorem) one checks 
immediately that 

= T7 1 1 + ^4^-IT A + £ 2^n(C-)] v + 0(e 2 ) 

7 + L 3d 7 + 7 + 

(the nonlocal operator II does not appear in the final equations because of 
the identity V • UV = V • V for all V G H 1 (R d ) d ). 

Step 2. The case a± > 0, /3 = 0, o<2 = 0. To use the BBM-trick, remark 
that for all a± > 0, 

V • v = (1 - ai)V • v - 01(7 + 8)d t ( + 0(e). 

Substitute this relation into the third-derivative term of the first equation 
of the system derived in Step 1. 

Step 3. The case a\ > 0, (3 > 0, 02 = 0. It suffices to replace v by 
(1 — ^/3A)v j g in the system of equations derived in Step 2. 
Step 4. The case a± > 0, (3 > 0, a.2 < 1- This is exactly as in Step 4 of 
Theorem El □ 

3.2 The Shallow Water/Shallow Water Regime: j»~/j 2 <l 

Contrary to the regimes investigated above, large amplitude interfacial de- 
formations are allowed for both fluids, as e ~ £2 = 0(1). As in the previous 
section, an asymptotic model can be derived from (|14p by replacing the 
operators G^feC] and H M ' 5 [eC] by their asymptotic expansions, provided by 
Proposition [2] and Corollary H] in the present regime. The following theorem 
shows that the internal wave equations are consistent in this regime with 
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the Shallow water /Shallow water system, 

' d t C + • (fciQ[^eJC](fov)) = 0, 

< 5 t v + (l-7)VC 

+ §v(|v - ^0[^^C](^2V)| 2 - (^IHI^CK^v)! 2 ) = 0, 

(29) 

where hi = 1 — e£, h% = 1 + e5£, and the operator is defined in Lemma 

Theorem 4. Let < <5 min < 5 max < (1 - 5(1 - Hi)) -1 . The internal waves 
equations p4\ ) are consistent with the SW/SW equations [2§\) in the sense of 
Definition^ with a precision O(fi), and uniformly with respect to e G [0, 1], 
A* G (0,1) and 5 G [<5 min ,5 ma!B ]. 

Remark 26. Taking 7 = and 5 = 1 in the SW/SW equations (f29|) yields 
the usual shallow water equations for surface water waves (recall that it 
follows from Lemma [3] that V- [(l-eC)0[-eC]((l + eC)v)] = V- ((l + eC)v)). 

Remark 27. In the one-dimensional case d = 1, one has 

— TT^^^CK^v) = ^ h2 - 
7 + <5 7 + 5 5/11+7/12 

and the equations (|29p take the simpler form 

' ftC + ft,(»gftfev)=o, 

' ^v + (l- 7 )^C + |54g^g||vP)=0, 

which coincides of course with the system (5.26) of [T7J. The presence of 
the nonlocal operator £2, which does not seem to have been noticed before, 
appears to be a purely two dimensional effect. 

Proof. First remark that with the range of parameters considered in the 
theorem, one has [i ~ [12 as [i — > while e ~ ei = 0(1). 
By the definition (|16p of v and using Proposition [2] and Corollary HI one 
deduces from (fT4"|) that 

/ ^C-V-((l-eC)VVi) = 0(/i), rsm 
\ diy + (i _ 7 )VC + f V(|H^[ £ C]^i| 2 - 7 |V^i| 2 ) = 0(/x). ^ 

Recall now that H^' 5 [e£]^i = v + 7V Va! since moreover one also gets from 
Corollarydthat H. 11 ' 6 [e(]ipi = -5£2[e 2 C](/iiV'0i)+O(/x), it is straightforward 
to deduce that 

v + 7 vv>i = -<5n[£ 2 cPiW>i) + o(^) 
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Multiplying this relation by h 2 and taking the divergence, one gets 

V • (h 2 v) + 7 V • (/12VV1) = -<JV-(M[£2C](fciV^i)) + V-0(M) 

= -5V-(/iiVVi)+V-0(/x), 

where the second equality comes from the definition of the operator Ofs^C]- 
We thus have 

v • (C 1 + ^^OVVO = — • (h 2 v) + v • OM, 
7+0 7+0 

and we can therefore use Lemma [3] to conclude that 

7+0 7+0 

and consequently 

= v + 7 V^i 

= V -^T7 £J [^T7 £ 2C](^V)+0( M ). 

7+0 7+0 

Replacing V^i and H^' 6 [e£]ipi by these two expressions in (f30|) yields the 
result. □ 



3.3 The Shallow Water/Small Amplitude Regime: fi <C 1, 

£ 2 < 1 

Derived here are various models corresponding to the case when the upper 
fluid layer is shallow, but this restriction is not required of the lower layer. 
The interfacial deviations are thus not necessarily small relative to the the 
upper fluid depth, but they are small relative to the undistrubed depth of 
the lower layer. Different systems of equations obtain, depending on the 
sizes of the paramaters e, // and S (and thus e 2 and fi 2 ). 

3.3.1 The Shallow Water/Full Dispersion Regime: /U ~ e| <g; 1, 

e ~ \i 2 ~ 1 

In this regime, the internal waves equations are consistent with the Shallow 
Water /Full Dispersion system, 

dt( + -V • (fciv) - • (hi\D\ coth(^|D|)n(/i lV )) = 0, 

J 7 7 

d t v + (1 - 7 )vc - y v V I 2 " 2 ^ v • coth(v^|D|)n(/nv))] = 0, 

(31) 
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where h\ = 1 — eC, and II 



vv 2 



Theorem 5. Let < c mm < c max and fi% im < [12 < 1^2 ax ■ The internal 
waves equations fi!4\ ) are consistent with the SW/FD equations $31\) in the 
sense of Definition^ with a precision O(fi), and uniformly with respect to 
e G [0,1], /i£ (0, 1) and 5 G (0, 1) satisfying the conditions 

c min <4^< c max and iit n < 4> < V2 ax - 
e z o z o A 

Remark 28. The SW/FD system (j31j) . which as far as we know is new, is a 
generalization of the results of §5.4 of [T7] to the two-dimensional case d = 2 
and to the case of a lower layer of finite depth (the case of an infinite lower 
layer is formally recovered here by taking T w = 1 in (|31D ). 

Proof. First remark that with the range of parameters considered in the 
theorem, one has e\ ~ \i and e ~ \ii ~ 1 as \i — > 0. 

Proposition [2] implies that ^G^eQ^i = V • (hiVipi) +0(ji) while it follows 
from the definition of v and Corollary [5] that 



1 

— v -\ - ■ 

7 7 2 T M2 

One then concludes the proof exactly as in the previous sections. □ 
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3.3.2 The Intermediate Long Wave Regime: /i. ~ e ~ e% -C 1, 

^2 ~ 1 

In this regime, a one-parameter family of intermediate long wave systems 
may be derived from the internal waves equations. These depend upon the 
parameter a and have the form 

( ct 

[1 + V/HD coth(V££|I>| MC + -V • ((1 - eC)v 
7 __ 7 

-(l-a)^|L>|coth( v ^| J D|)V-v = 0, (32) 
5 t v+(l- 7 )VC-^V|v| 2 = 0. 

Theorem 6. Let < c min < c max , fif in < fi 2 < nf ax - The internal 
wave equations fll^| ) are consistent with the ILW system \3%fy in the sense of 
Definition^ with a precision O(p), and uniformly with respect to e G [0, 1], 
fi G (0, 1) and 5 G (0, 1) satisfying the conditions 

c min < _M < c max Qnd < M < ^max_ 
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Remark 29. In dimension d = 1 and with a = 0, (|32p corresponds to (5.47) 
of [T7] . However this system is not linearly well-posed. It is straightforward 
to ascertain that the condition a > 1 insures that f|32|) is linearly well-posed 
for either d = 1 or c£ = 2. 

Remark 30. The ILW equation derived in |23[ [25] is obtained as the uni- 
directional limit of the one dimensional (d = 1) version of (|32|) - see for 
instance §5.5 of [T7] . 

Proof. Step 1. The case a = 0. We are working with the regime fj, ~ e 2 ~ 
£2 <C 1 and H2 ~ 1 as — > 0. In this situation, Proposition [2] allows us 
to write ^G^eQ^x = V • ((1 - eC)V</>i) + 0(/i) while it follows from the 
definition of v and (1241) that 



VVi = --v + 4^v + 0(m). 
7 T T w 

Substituting these two relations into the internal wave equations (|14p leads 
to the advertised result with a = 0. 

Step 2. The case a > 0. This result follows from Step 1 and the observation 
that 

V • v = (1 - q)V • v - a-fd t C + 0(e, y/JT). 

As mentioned already, the restriction on a is not to obtain consistency, but 
rather to ensure linear well-posedness. □ 

3.3.3 The Benjamin-Ono Regime: [i ~ e 2 <S 1, fi2 = 

For completeness, we investigate the Benjamin-Ono regime, characterized by 
the asumption 5 = (the lower layer is of infinite depth) . Taking /12 = 00 in 
(|3"2j) leads one to replace coth(y / 7i2"|rJ|) by 1. The following two-dimensional 
generalization of the system (5.31) in [T7j emerges in this situation. 

f [1 + V^f \D\]d t C + ±V • ((1 - eC)v) - (1 - a)^\D\V • v = 0, 
1 5 t v+(l- 7 )VC-^V|v| 2 = 0. 

Neglecting the 0(yfii) = 0(e) terms, one finds that £ must solve a wave 
equation (with speed ■J^~^)- Thus, in the case of horizontal dimension 
d = 1, any interfacial perturbation splits up at first approximation into two 
counter-propagating waves. If one includes the 0(y / 7i, e) terms, one obtains 
the one-parameter family 

(1 + sfJi-\d x \)d t C, + cd x ( - e-cd x C 2 - ^c(l - 2a)\d x \d x ( = 0, (34) 
7 4 27 
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of regularized Benjamin-Ono equations (see ??). Here, c = \J~^~- The 
usual Benjamin-Ono equation is recovered by taking a = 0. 

A Proof of Proposition [3] 

The proof is made in five steps. 

Step 1. Coercivity of the operator V^ 2 z ■ Q^[s 2 (}V£ z . Exactly as in Prop. 
2.3 of [2], one may check that 

veei^ 1 , e-Q w [e 2 C]© > tI@I 2 5 

with k = fc(^,e v ^t,e 2 |Clw 1 .°°) > °- 

Step 2. Existence of a unique solution to jfgjp. Owing to Step 1, existence of 
a solution and uniqueness up to a constant is provided by classical theorems 
(e.g. Section V.7 of [34]). provided that the source terms and Neumann 
conditions satisfy the compatibility condition 

/ ^xz ■ h = f ( v^ v ■V + e.-h)- f e z • h. 

' ./{z=0} J{ 2 =-1} 

This latter restriction is valid in the present circumstances on account of the 
divergence theorem. 

Step 3. L 2 -estimate on z u. Multiplying (i2~Tj) by u, integrating by parts 
on both sides, and using the Neumann conditions leads to 

/ V£> • M V£> = - / V ■ ^Nu + / h • V£>. 

JS J{z=0} Js 

A direct consequence of Step 1 and the Cauchy-Schwarz inequality is the 
inequality 

||V£>|| 2 < fe(||h|| \\V%u\\ + \V\ H ^\^u\ H -y,). 

It follows from the trace theorem that 

\^Vu\ H - 1/ 2 < Cst (||V/^Vtt|| + HA^v^&Vull) 
< Cst(||V^>|| + v^||V^>||). 

It is concluded that 

l|V£>|| < C(±-,eT\^ ax M\w^){\U + \V\ H m)- 
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Step 4. H s -estimate (s > 0) on V^ z u. Let v = A s u. Multiplying §21$) by 
A s on both sides, it results that v solves the system 

| V^-Q^[ £2 C]V^ = V^-h, in S, 

\ d n v u=0 = ■ A S V + e z • h\ z=0 , d n v u ^_ L = e z • h u= _ 15 



with h = A s h + [Q^ [e 2 (], A s ]V^ z u. From Step 3 and the definition of v, it 
is thus deduced that HA^Vj^iiH is bounded from above by 

C(-^,er x ,»r x , |ClH/i.~)(HA s h|| + \V\ HS+ y, + ||[Q W M, A S ]V£>||). 
Using the expression for Q^^C] an d the commutator estimate 
\[A s ,f]g\ 2 <C\Vf\ 

which holds for some constant C which depends upon s > — % and to > f 
(see Th. 6 of [23); we obtain 

We thus get an estimate on ||A s V^?.,u|| in terms of ||A s_1 V^? z u|| which, 
together with Step 3 (i.e. s = 0) allows us to derive the following relation 
by induction (and interpolation when s G (0, 1)): 

v s > o, \\a s v% z u\\ < c(^, e r x ^rMCI^oo + ^ + i})(l|A s h||+|y| HS+1/2 )- 

Step 5. H s ''-estimate (s > 0) on c^V^u. First remark that using the 
equation yields 

1 + ^ U + 1) W 

1 + e 2 C ' 
-VJ^V • ((1 + e 2 C)^Vu - Jjle{z + l)V(d z u) 

wc\ 2 

+v^/iVC • (V^Vn) - 2fie 2 (z + l)±-^-d z u, 

t + £2<, 

from which one obtains the estimate 

||A S 5^|| < C( £ ^,//^,|C|^ max Oo +2 . s+ i>)(l|A s V^h|| + V^I|A s+1 V^^||). 
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Use this to write 

||A s cW£>|| < ^\\A s d z Vu\\ + \\A s d 2 z u\\ 



< c^^rMci^ oo +2 , s+ i})(iiA s v^ i2 h|| +v^iia s+i v^ 

With the help of Step 4, one obtains the inequality 

\\A s d z v% z u\\ < c(i, £ r s ^rMci^ {t o +2 , s+ i})(iihii^ +1 , 1 + \v\ HS+m ). 

Step 6. Conclusion. By the trace theorem we may assert that for all s > 0, 

Cst 

|V«i |/f S < Cst \\Vu\\ H s+i/2,i < ^=||V^%n|| Hs+ i/ 2 ,i- 



The desired result now follows from Steps 4 and 5. 
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